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PREDICTION OE PROPERTIES OE I NT RA PLY HYBRID COMPOSITES 


by C. C. Chamis* and J. H. Sinclair* 


National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 


original PAGE IS 
POOR QUALITY 


Equations based on the mixtures rule are presented for predicting the 
physical, thermal, hygral, and mechanical properties of unidirectional intraply 
hybrid composites (l H1C) ffoin the corresponding properties of their constituent 
composites. Bounds were derived for uniaxial longitudinal strengths, tension, 
compression, and flexure of U1HC. The equations predict shear and flexural 
properties which agree with experimental data from UIHC. Use of these equa- 
tions in a composites mechanics computer code predicted flexural moduli which 
agree with experimental data from various intraply hybrid angleplied laminates 
(RIAL). It is indicated, briefly, how these equations can be used in conjunction 
with composite mechanics and structural analysis during the analysis/design 
process. 


INTRODUCTION 

The significance of fiber hybrid composites as an emerging structural 
material is extensively discussed in a recent state-of-the-art review (ref. 1). 
The various types of hybrid composites that have been used and their special 
advantages and disadvantages are covered in considerable detail in that review. 
The available methodology for analysis and design as well as areas that need 
further research are also covered. A need for equations for predicting the 
various properties of intraply hybrid composites based on the corresponding 
properties of their constituent composites was pointed out. The properties of 
interst include the following: physical, thermal, hygral, and mechanical . The 
objective of this paper is to present approximate equations which are based on 
the mixtures rule and also satisfy the conditions of mechanics for predicting all 
these properties. These equations ( then be used in conjunction with com- 
posite mechanics computer codes for designing structural components made 
from intraply hybrid composites. 

Aerospace Engineers, Material and Structures Division, NASA Lewis 
Research Center, Cleveland, Ohio 44135. 
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The model and assumptions which justify the use of the method as well as 
the equations for predicting the various properties are described in detail. 
Comparisons with limited available experimental data are made to assess the 
adequacy ol the method. I*, is indicated how these equations can be used in con- 
junction with available composite mechanics and structural analysis computer 
codes tor the ana lysis/de sign of structural components made from intraply hy- 
brid composites. The symbols used In the equations are defined when they ap- 
pear and they are also summarized in appendix A for convenience. 


MODEL AND ASSUMPTIONS 

The volume average (rule of mixtures) method for predicting approximately 
the various properties of unidirectional intraply hybrid composites (UIHC) Is 
based on the model of a UIHC depicted schematically in figure 1. In this figure 
the model is represented by an array of parallclpiped cells which are differenti- 
ated by either empty or filled circles. The empty circles represent primary' 
(majority) fibers or fILcr bundles while the filled circles represent secondary 
(hybridizing) fibers or fiber bundles. Rectangular coordinate system axes 
(1, 2, 3) are defined where the 1-axis is parallel to the fiber direction, 2 is 
along the width (transverse direction! and 3 is through the thickness. 

The assumptions required to apply the method of volume averages for pre- 
dicting the various properties of UIHC are as follows: 

(1) Each pai allelpiped cell is homogeneous and its properties are the same 
as those of either the primary fiber composite or the secondary fiber composite. 

(2) The cell boundaries are contiguous and there is a complete bond between 
cells. 

(3) The various properties of the constituent composites are known either by 
experiment or by use of composite micromechanics. 

(4) The UIHC and its constituent composites exhibit linear behavior in the 
regions of interest for this investigation. 

(5) Plane sections remain plane after deformation resulting from the applica- 
tion of uniform heat flux, moist environment, or stress on planes parallel to the 
coordinate planes (fig. 1). 

A direct consequence of assumption (5) above is that the cells in the model 
(fig. 1) are connected in parallel in all three coordinate axes planes. For exam- 
ple, the stress in each cell due to a uniform stress in the 2-direction (transverse 
direction) is proportional to its stiffness. It is generally accepted that the rule 
oi mixtures applies to properties along the fiber direction because along this 
direction the elements in the model are connected in parallel. However, the ap- 
plication of the rule of mixtures in the transverse and through the thickness 
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directions is not so obvious. As can be seen from figure 1, the elements in 
these directions are connected in both parallel and series. Another consequence 
from the above assumptions is that the I'lllC is homogeneous as viewed from the 
outside. In view of the ubovc discussion, the properties of UIHC can be approxi- 
mately given by .he rule of mixtures which can be expressed by the following 
general equation ; 


P HC _ P PC + V SC (P SC " P PC ) 


( 1 ) 


where P denotes the property of interest, V denotes hybridizing (secondary) 
composite volume ratio, subscript 1IC denotes hybrid composite property, 
subscript PC denotes primary composite property, and subscript SC denotes 
secondary composite property. Numerical subscripts denote direction. Exten- 
sion of the equation for intraply hybrids from more than two different fibers is 
achieved by adding terms analogous to the second term in the right hand side of 
the equation (e.g., using superscripts the second term would be Vg^(Pg^ - Pp^,), 
etc . ) . 

Equation (1) expresses each property of the U1HC as a linear combination of 
the corresponding properties of the constituent composites and satisfies both 
boundary conditions at Vg q equal "0" and 1 . This equation, when applied to 
mechanical properties, atisfies the three conditions of mechanics (force equilib- 
rium, compatibility, and stress-strain relationships). Typical properties of 
various unidirectional fiber composites which have been used to fabricate UIHC 
are summarized in Table I. Unidirectional composite properties which are not 
readily available and are usually predicted by using composite micromechanics 
are summarized in Table II. 


EQUATIONS FOR THE VARIOUS PROPERTIES 

The properties of UIHC that may be required for structural and stress analy- 
ses of UIHC include physical, thermal, hygral, clastic, and strength and are sum- 
marized in appendix B. Selected equations are described below to clarify the no- 
tation used and to illustrate significant points. The desired equations are obtained 
by appropriate substitutions in equation (1). Constituent composite properties for 
use in the equations may be obtained from Tables I and II or can be determined 
from fiber and matrix properties using composite micromechanics (refs. 2 to 4). 
The various equations that may be needed for structural and stress analyses due 
to thermal hygral and mechanical loads are summarized in appendix B. 
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Thermal, llygral, and Elastic Properties 
The heat conductivities oi a UIHC are given by 

K IlCi = K PCi + V SC (K SCi " K PCi ) 

where the subscript l is a generic index and denotes direction 1, 2, or 3 along 
which the heat conductivity is desired. The thermal expansion coefficients, 
moisture diffusivities, moisture expansion coefficients, normal moduli, and 
flexural moduli of UIHC are given by similar equations where K is replaced by 
the appropriate symbol for the desired property. The shear moduli of UIHC are 
given by 

G IICiJ = G PCij + V SC^ G SCij ' G PCiJ ) (3) 

where the subscripts l and j denote planes 1-2, 2-3, or 1-3 in figure 1. The 
corresponding Poisson's ratios are given by 

^HCij ‘'PCij + V SC^SCiJ “ ‘'pCiJ* (4) 

where the subscripts ij denote contraction in the j direction due to stress in 
the i-direction. 

Unidirectional fiber composites (thoroughbreds) are generally assumed to be 
transversely isotropic, that is, the properties in the 3-direction are equal to cor- 
responding properties in the 2-direction (fig. 1). For example: E^ - E 2 ; - 

Gj^* an< * = Properties associated with the 1-3 plane (G 23 and v 2 ^) 

are not readily measured. Experimental value for these properties are virtually 
nonexistent. These properties are usually estimated by using composite micro- 
mechanics (ref. 2). As was already mentioned, the properties summarized in 
Table 11 were predicted by using the composite micromechanics available in the 
computer code (ref. 4). 


Strengths 

The strengths (fracture stresses) (S) of interest in performing stress analy- 
ses are: longitudinal tensile and compressive, transverse tensile and compres- 
sive, intralaminar (in-plane), interlaminar (short beam) and through-the- 
thickness shear; and longitudinal and transverse flexural (bending). The equa- 
tions for predicting approximately the strengths of UIHC are obtained by appro- 
priate substitutions in equation (1) as was the case for the properties already 
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described. Typical strengths which are readily measured, for constituent uni- 
directional composites are summarized in Table 1. Those predicted by using 
composite micromechanics ure summarized in Table 11. 

The longitudinal tensile strength of a L'lHC is given by (refer tc fig. 1) 

S HC1T S PC1T + V SC (S SCIT “ S PC1T ) (5) 

where I denotes strength along the fiber direction and T denotes tension. 

The other normal strengths are obtained by substituting the appropriate sub- 
scripts in equation (5). 

The intralaminar shear strengths are given by 

S HC12 = S PC12 + V SC (S SC12 “ S PC12 ) 

S x»C23 S PC23 + V SC (S SC23 * S PC23 ) (7) 

where the subscripts 12 and 23 denote the shear planes (fig. 1). The inter- 
laminar shear strength (short -beam-shear) of UIHC is given by 

S 11CSB S PCSB + V SC (S SCSB ” S PCSB ) (8) 

where the subscripts SB denote short-beam shear. The longitudinal flexural 
strength is given by 

S HC1F S PC1F + V SC (S SC1F " S PC1F ) (9) 

where the subscript F denotes flexural property. The transverse flexural 
strength is obtained by replacing subscript 1 by 2 in equation (9). 

One important point need be made in connection with equations (5) to (9). 

As was already mentioned the equations show that the strength of the UIHC is a 
linear combination of the strengths of its constituent composites. Intraply hy- 
brid composites which are of good quality, thoroughly mixed and completely 
bonded, should exhibit strengths which fall on or above the straight line de- 
scribed by the corresponding equation (full hybrid response). 

Occasionally points will fall below the line. The causes for this are prob- 
ably poor quality fabrication and/or a bad batch of material which results in 
ineffective use of the intraply hybrid. 
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Lower Bounds on Longitudinal Strengths 


Equations (5) and (9) will overestimate the longitudinal tensile 
compressive (Sj 1C1 ^) and flexural ( s jj^; j. j.') strengths in UIHC with partial hy- 
brid response (incomplete bond). Equations for predicting the lower bound 
strengths in such IIHC can be derived by assuming that the primary composite 
induces fracture at low values of the hybridizing ratio, V SC* and that the sec- 
ondary (hybridizing) composite induces fracture at high values of Vg£. The 
resulting equations are: 

S HCly c E HC1 for V SC s V T (10) 

^PCl 


c c V S 

a HCly N SC a SCly 


for V SC “ V T 


( 11 ) 


where 



1 

L | S SC1~) £ SC1 

S PCly E PC1 


( 12 ) 


where y denotes T, C, or F lor tension, compression or flexure, respec- 
tively, and Vj denotes the hybridization ratio (volume ratio of secondary com- 
posite) at which transition of the fracture mode occurs. Implicit in the deriva- 
tion of equations (10) to (12) is the assumption that the fracture strain of the 
libers in the primary composite is smaller than that of the fibers in the second- 
ary composite. If the reverse is true, then the role of the two composites in the 
equations should be interchanged. The graphical representation of equations (5), 
(9), (10), and (11) is shown in figure 2 for a graphite fiber/glass fiber intraply 
hybrid, HMS/S -G(901-S)(data from Table I). Partial hybrid response (incom- 
plete bond or nonuniform mix) results in ineffective utilization of the intraply 
hyfcnci. 


Other Properties 

Equations for predicting IIHC properties such as impact resistance, fatigue, 
and fracture toughness, can be derived by appropriate substitutions in equa- 
tion (1). These equations, however, would only be applicable to initial damage 
wnere the intraply hybrid composite responds like a homogeneous material. 
Consequently, these equations will not account for cases in which a combination 
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of failure modes participates to propagate flaws. Therefore, equations for im- 
pact, fatigue, and fracture toughness are not given herein. 


COM PA RISC NS 

Only limited experimental property data are available where the same com- 
posite system was used for both constituent composites and I'lHC. Comparable 
data for thermal and hygral properties are not available. Comparisons between 
predictions and available experimental data are summarized in Table 111 for 
interlaminar shear strength (short beam shear), longitudinal flexural strength, 
and flexural modulus for three different intraply hybrids. The experimental 
data for these comparisons are from reference 5. The data for the relevant 
constituent composites is summarized in the top part of the table. The predic- 
tions were obtained bj using equations (B8), (1UH), and (B20). As can be seen 
the predicted values are within 12 percent of the interlaminar shear strength, 

17 percent of the longitudinal flexural strength and 11 percent of the longitudinal 
flexural modulus. These results are considered very good for the following 
reasons: 

1. The interlaminar shear and flexural responses are complex structural 
responses since they include combinations of tensile, compressive and shear 
stresses and strains. 

2. The experimental data are evenly distributed above and below the theo- 
rectical predictions which is consistent with the observations made earlier with 
respect to complete and incomplete bonding. 

3. The types of intraply hybrids investigated in reference 5 are the "tow-by- 
tow" variety. In these types of hybrids the hybridizing (secondary) fibers are 
not as uniformly dispersed as was assumed in the physical model. 

The comparisons are also presented graphically in figure 3 together with 
the bounds for full and partial hybrid response. Note the interchange of the role 
between the primary and secondary composites for the AS/KEV L1HC. As can 
be seen in figure 3 the experimental points fall either between the bounds or 
above the full hybrid response bound. The largest percent difference between 
predicted and measured values (17 percent for the longitudinal flexural strength 
of the HMS/S-G intraply hybrid) is well above the partial hybrid response range 
which is about 30 percent relative to full hybrid response (fig. 3). 

The important conclusion from these comparisons is that the rule-of- 
mixturcs equations predict intraply hybrid mechanical properties which are in 
good agreement with measured data for hybridization volume percent up to about 
30. An interesting observation from the data in Table III is that the intraply 
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hybrid AS/KEV shows an apparent synergistic effect since its flexural strength 
oi 2u:i ksi is about li percent higher than that of its constituent composites, which 
is 275 ksi. 


APPLICATIONS TO DESIGN 

The design ol structural components made from anglcplied intraply hybrid 
composites requires analyses which use composite macromechanics, laminate 
theory, and suitable failure criteria. I sc of composite macromechanics, lami- 
nate theory and fuilurc criteria involve lengthy tedious computations. These are 
generally done via composite mechanics computer codes. One such code is de- 
scribed in reference -1 us was already mentioned. The composite mechanics 
codes are used twice during one structural analy sis and/or design iteration cy- 
cle. In the first step, they are used to generate the properties required for 
structural analysis such as axial and bending stiffnesses. In the second rep, 
they are used to calculate the ply stresses and assess the luminatc strength 
using the loads and/or the displacements determined from the structural analy- 
sis. Structural analy ses to determine design variables such as displacements, 
forces, vibrations, buckling loads, and dy namic responses include application 
of corresponding special areas ol structural mechanics lor simple structural 
elements. General purpose finite element programs such as NASTRAN are used 
for the structural analysis of complex structural shapes, large structures made 
from simple structural elements, and structural parts made from combinations 
of simple elements such as bar, rods, and plates. 

Composite mechanics codes require unidirectional composite (ply) proper- 
ties as input. These properties for 1 1 1 1C can be determined from the equations 
described previously . Once these properties are available the structural analy- 
sis and/or design of structural components made f ram intraply hybrids is the 
same as lor "thoroughbred” composites. As an example, the equations were 
used to provide the required input to a computer code (ref. 4) to predict the flex- 
ural (bending) modulus lor several intraply hybrid anglcplicd laminates (111AL). 

The results obtained are compared to the experimental data in Table 1V ? . As can 
be seen the predicted values are in reasonably good agreement with the measured 
data. This reasonably good agreement for the flexural moduli of 1HAL is con- 
sidered to be a further validation of the rule of mixtures method since the input 
to the code required four properties a '>d •'em") P ru " 

dieted by four different equations lor each I'illC. Aside from the validation as- 
pect flexural moduli were calculated because they are important in determining 
structural responses (behavior variables) such as bending displacements, buckling 



loads, vibration frequencies, and local deformations due to out-of-plane point or 
impact loads. 

The previous discussion leads to the following conclusion and anticipation. 
The equations described herein predict mechanical properties at both the uni- 
directional composite and laminate levels which are reasonably good agreement 
with experimental data from three different intraply hybrid composite systems, 
because of this good agreement, it is antici|>ated that these equations will apply 
equally well to other properties for these intraply hybrids and to properties for 
Intraply hybrids from different constituent composite systems. 


SUMMARY 

Equations wure presented, based on the rule of mixtur* s, for predicting 
physical, thermal, hygral, and mechanical properties (including strength) of 
unidirectional intraply hybrid composites (L'lHC). Results predicted by using 
these equations are in good agreement with lbnited available data for L'lHC. 
Experimental strength data fell within or above bounds predicted by equations 
derived herein for this purpose. Use of these equations with a composite me- 
chanics computer code, predicted flexural moduli of intraply anglcplicd lami- 
nates (1HAL) which were in reasonably good agreement with experimental data. 
The availability of these equations makes it possible to design, analyze, and 
optimize structural components from IHAL based on a large variety of available 
constituent composites . 


APPENDIX A - SYMBOLS 

C heat capacity 

D moisture diffusivity 

E normal (elastic) modulus 

G shear modulus 

IHAL intraply hybrid anglcplicd laminate 

K heat conductivity 

P property - physical, thermal, hygra! or mechanical 

S strength 

L'lHC unidirectional intraply hybrid composite 

V volume (hybridization) ratio (volume of hybridizing (secondary) 

composite /I otr! volume of intraply hybrid composite) 
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V,j. volume ratio at which transition of fracture modes occurs (eq. (12)) 

1,2,3 reference coordinate axis (fig. 1) 

u thermal expansion coefficient 

(i moisture expansion coefficient 

v Poisson's ratio 

p density 

Subscripts: 

C composite property, compression 

F fiber property, flexure 

H hybrid composite property 

i generic index (1, 2, or L) 

j generic index (1, 2, or 3) 

P primary composite 

S secondary composite 

SB short beam (interlaminar) shear 

T tension, transition 

1,2,3 direction associated with corresponding coordinate axis 
y generic index (T, C, or F) 

APPENDIX B - UST OF EQUATIONS 

Several of the equations which are obtained by using the rule of mixtures 
(appropriate substitutions in eq. (1) in the text) and which may be used to predict 
various properties of intraply hybrid composites are summarized below. The 
notation is defined in appendix A. 

Physical properties - fiber volume ratio and density: 

V HCF = V PCF + V SC (V SCF ~ V PCF ) (B1) 


» 


P HC ~ P PC + V SC (p SC " P PC ) 


(B2) 
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Thermal properties - heat capacity, heat conductivities, and expansion coeffi- 
cients: 


C HC C PC + V SC (C SP _ C PC ) 

(B3) 

HCi - K PCi + V SC (K SCi " K PCi ) 

(B4) 

MCI a PCi * V SC (a SCi " Q PCi ) 

(B5) 


where i is a generic index and denotes direction 1, 2, or 3 along which the 
property is desired (fig. 1). 

Ilygral properties - moisture diffusivlties and moisture expansion coefficients: 


U HCl D PCi + V SC (D SCi " 

U PCi ) 

(BO) 

^HCi * ^PCi + V SC ( ^Cl " 


(B7) 

Elastic properties - normal and flexural moduli, 
ratio : 

shear moduli and Poisson's 

E HCi " E PCi + V SC (E SCi “ E PCi ) 

i- 1.2,3 

(BO) 

G HCij " G PCij + V SC (G SCij " G PCij ) 

ij * 12,23,13 

<B9) 

^llCij = ^PClj + "sC ( *'PCij " *'PCij ) 

ij = 12,23,13 

(BIO) 

E HCiF 3 E PC1F + V SC< E SCiF “ ^PCiF 1 1-1,2 

(Bll) 


Strengths - longitudinal tension and compression, and transverse tension and 
compression: 


S HC1T = S PC1T + V SC (S SC1T " S PC1T ) 
S HC1C = S PC1C + V SC (S SC1C " S PC1C ) 
S HC2T S PC2T + V SC (S SC2T " S PC2T* 


(B12) 

(B13) 

(B14) 


S 1IC2C ~ S PC2C + V SC (S SC2C “ S SC2C ) 


(B15) 
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Strengths - intralaminar (in-plane) shear and through-the-thickness shear: 

S HC12 = S PC12 + V SC (S SC12 " S PC12> 

S HC23 “ S PC23 + V SC (S SC23 " S PC23 ) 

Strengths - flexural and interlaminar shear (short-beam-shear): 

S HC1F * S PC1F + V SC (S SC1F ' S PC1F ) 

S HC2F = S PC2F + V SC (S SC2F " S PC2F ) 

S HCSB = S PCSB + V SC (S SCSB ’ S PCSB ) 
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TABLE U - UNIDIRECTIONAL COMPOSITE PROPERTIES PREDICTED BY COMPOSITE MICROMECHAXICS 


USING THE COMPUTER CODE REF 4 



— 

Boron 

Mod-1 

Mod-U 

T3M 

T75 

Kevlar-49 

<KEV» 

S-glass 

tS-G» 

E -glass 
|E-C| 

Heat capacity, Btu/(lb/ > F) 

0.29 

0 20 

0.20 

0.2 

0.2 

0.21 

0.19 

0.19 

Heat co nductivity (Kj>, Btu/<hr)(ft 2 )(°F/in) 

12.0 

291 

291 

348 

291 

4.6 

4.6 

4.5 

Heat conductivity (K.„ K^), B»u/<hr><ft 2 )(°F/»n) 

4.0 

3.7 

3.7 

5.2 

3.7 

2.7 

2.7 

2.5 

Thermal expansion coefficient (o^), 10~ 6 in/in/ > F 

3.1 

-0.2 

0.01 

0.5 

-0.2 

-1.4 

3.9 

4.1 

Thermal expansion coefficient (to,, a 3 ), 10 -(i in/in/ > F 

16.3 

23.2 

23.2 

16.2 

23.4 

32.2 

16.1 

16.0 

Shear modulus (G 2a ). 10 psi 

0.58 

0.30 

0.33 

0.37 

0.34 

0.24 

0.50 

0.49 

Poisson's ratio (iv 

0.53 

0.40 

0.41 

0.47 

0.41 

0.42 

0.50 

0.49 

Through-the-tliickness shear strength <S 23 ), ksi 

13.2 

6.0 

7.7 

8.5 

8.2 

8.3 

8.2 

6.1 


Note: Subscript l denotes property along the liber direction. 

Subscript 2 denotes property transverse to the fiber direction 
Subscript 3 denotes property through the thickness. 

Double subscripts denote properties in that plane, see also fig. I. 
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TABLE IU - COMPARISONS OF EXPERIMENTAL AND PREDICT Ft) PROPERTIES 


(Experimental data Irom ref. 5.| 


Constituent composite* 

Rho rt -beam - (thru r 
strength, k»i 

Flexural strength, kai 

Flexural modulus, 

10° pal 

Experi- 

mental 

Pre- 
diet t^d 

Differ- 
ence, % 

Experi- 

mental 

Pre- 

dicted 

Difler- 
ence, % 

Experi- 

mental 

Pn- 

dieted 1 * 

Dlfter- 
enoe, *7 


AS 

hi 

g 

g 

275 



18 




I1MS 

WSm 



172 



28 




S-O 




275 



» 




KEV 

mm 

B 

B 

98 


B 

11 

B 

B 

lnlrupl) hybrids 





g 





Hybrid 

Con stituents, 'I 

■ 




g 





HMS/S-G 

73.5/26.5 

n 

9 

12 

171 


17 

21 

23 

10 

ASS-0 

75. 0/25. 0 

EM 

19 

0 

El 

131 

1 

18 

16 

-11 

AS/ KEV 

67 .5/12.5 


16 

7 

MM 

ESI 

•14 

18 

17 

-6 


“Predictions bused on full hybrid rvsponae. 


TABLE IV. - I NT HA PLY HYBRID COMPOSITE FLEXURAL MODl'LUS (10° pal) 
COMPARISON BETWEEN MEASL RED AND PREDICTED DATA 


Ply configuration 

AS/S-G 

AS/ KEV 

HMS/S-G 

Measured 

Predicted 

Measured 

Predicted 

Measured 

Predicted 

( >40,0, 10,0, -10) s 

7.75 

7.04 

7.95 

7.33 

6.81 

8.89 

( >22, 0, 22, 0, -22 )g 

12.3 

11.49 

12 35 

12.43 

12.71 

14.30 

( >45, 0, 45, 0, -45)g 

5.19 

5.29 

4.4 

5.38 

6.78 

6.71 

Fiber composition 

75/25 

90/10 

62/38 
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Figure 1. * Unidirectional intraply hybrid composite schematic. 
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Figure 2. - Effects ol partial hybrid response on the longitudinal 
strengths of unidirectional intraply h^trid composites 
iHMS/S-G iWl-Sll. 
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Figure }. Comparison ot predicted and experimental results lor llexural 
strength (experimental data Iron table III). 


